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Abstract 

The evolution of probability distribution functions (PDFs) of continuous density, velocity and 
velocity derivatives ( deformation tensor) fields in the theory of cosmological gravitational insta- 
bility are considered. We show that in the Newtonian theory the dynamical equations cannot be 
reduced to the closed set of Lagrangian equations. Since continuous fields from galaxy surveys 
need sufficiently large smoothing which exceeds the scale of nonlinearity, one can use the Zel'dovich 
approximation to describe the mildly non-linear matter evolution, which allows the closed set of 
Lagrangian equations. The closed kinetic equation for the joint PDF of cosmological continuous 
fields is derived in this approximation. The analytical theory of the cosmological PDFs with ar- 
bitrary (including Gaussian) initial statistics is developed, based on the solution on the kinetic 
equaiton. 

For Gaussian initial fluctuations, the PDFs are parametrized by only linear rms fluctuations a on 
given flltering scale. Density PDF P{p, t) and PDF M(Ai, A2, A3; t) of eigenvalues of the deforma- 
tion tensor fleld in the Eulerian space evolve very rapidly in non-linear regime. On the contrary, 
velocity PDF Q{v,t) remains invariant under non-linear evolution. For small a the Edgworth 
series is suggested to reconstruct P{p,t). Density PDF P{p,t) is close to the log-normal PDF for 
the CDM model at moderate a, but differs from that in general case. 



* Extended version of the Talk presented at the 9-th lAP Meeting "Cosmic velocity fleld", 
Paris, June 12-17, 1993. 



1 Introduction 



One of the standard metods to study the statistical property of the non-hnear Large Scale Struc- 
ture (LSS) of the universe is the two-point correlation function < y{fi)y{f2) >= ^i^K'^i,r2) and 



further, n-point functions ^, 



n, 



of the point-like statistical ensemble of galaxies, y may be 



the density, velocity components or other fields. The full statistics of the ensemble is described by 
the n-point PDFs P{yi, ■■■,yn) which are connected with ^("^ 



It is insufficient to know a few 



lower order cerrelation functions of galaxy distribution to define the full statistics of non-linear 
matter distribution in the universe. In this contribution we will consider the simplest one-point 
PDFs of the continuous non-linear cosmic fields. On the statistics of cosmological point processes 
see i. 

The theoretical development of the subject begun with the calculations of the lowest moments 
< 5^"'"' > of the density fields in the perturbation theory [0,[18|, in the Zel'dovich approximation 
|T^, for different cosmologies []TU[, for filtered field in numerical simulations for non-linear 
regime |]12| . 

Recently even more attention has been drawn to the PDFs which can probe the statistics of pri- 
mordial fluctuations. The standard CDM model with the scale free Gaussian initial fluctuations 
is hardly compatible with the set of observations. One way (fortunately not the only) is to switch 
the model to start with non-Gaussian primordial inhomogeneities. In the cosmic inflation scenario 
in principle, there is some room to introduce non-Gaussian fluctuations |2^, although it needs flne 
tuning in parameters, as a rule. Non-Gaussian perturbations also arise in other scenarios, where 
progenitors of inhomogeneities are topological defects, such as textures |^ , cosmic strings 
explosions 
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or 



The phenomenological attempts to design the density PDF at non-linear stage were made from dif- 
ferent arguments. In [JD| the formula for P{p) from "thermodynamical" treatment of gravitational 
system was suggested, which evokes some sceptizism. The log-normal mapping in order to mimick 
the non-linear density evolution was argued in |]l3l . Unfortunately the log-normal model does not 
work |T^ , although the features of the log- normal distribution have been recognized in the density 
PDF |2ll], [|l^, [^. An elegant analytical method has been developed ^ to summarize the series 
of moments for the small density fluctuations 5 in the approximation when < (5'-"^ >~ S'„(T^^"'~^\ 
In p6|, ||29| we calculated -P(p) and F{v) analytically in the "truncated" Zel'dovich approximation 
(i.e. with initial smoothing), and found a good agreement with the results of N-body simulations 
(with flnal smoothing). Recently two different methods from papers [0] and were compared 
resulting in a good mutual agreement [0. 

In |36] further approximation was made within the Zel'dovich theory to calculate P{p) after the fl- 
nal smoothing. Unfortunately the assumptions made are valid for small fraction of the Lagrangian 
volume and result does not flt the numerical calculations. The count-in-cell statistics which is 
tighly related to the P{p) was calculated numerically in [0]. In P{p) from a range of numerical 
simulations with different statistics was plotted. 

From observational side, the skewness measured from the IRAS QDOT |39] is in the range 
s/o" = 1 — 2, as expected from the Gaussian initial statistics. However, the IRAS galaxies are 
underrepresented in clusters, which correspond to high density tail of P{p), contributing to the 
skewness. This has little effect on P{p) itself because clusters occupy a small fraction of volume. 
The density and velociy PDF as deduced from 1.9 Jy IRAS survey and POTENT analysis is 
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consisient with Gaussian initial statistics [29|, the count-in-cell analysis of 1.2 Jy IRAS data also 
consistent with this hypothesis 0, although in all cases the larger sample volume is needed for 
more quantitative conclusion. 



2 Basic Non-Local Equations 



At the epoch of the LSS formation, most of the mass is in the form of dark matter of relic origin 
(for instance, the Cold Dark Matter), with no pressure; nonlinear LSS originated by gravitational 



instability of small initial fluctuations. Let x, v 



' dt ' 



p{t, x) and 0(t, x) be, respectively, the 



comoving coordinates, pecuhar velocity and density of dark matter (neglecting baryons), and 
peculiar Newtonian gravitational potential. In the Newtonian theory, the motion of dark matter 
before the particle orbits cross obeys a nonlinear system of equations 



at a 



0, 



(1) 



dv 
di 



a a 



(2) 

u 

Vj(l) = A7rGa\p-p), (3) 

where a{t) is a scalar factor of an expanding Universe, H = a/a and p is a mean density. These 
equations are valid in the single stream regime, and admit the evident generalization in the regions 
of multiple streams. 

For simplicity we assume that Inflation produces a flat universe = 1, and A = 0. The growing 
mode of adiabatic perturbations D{t) in the Einstein-de Sitter universe is D{t) = a(t) oc t"^^^. 
Let us perform further analysys of the basic equations. It is convenient to use the growing solution 
D[t) = a{t) as a new time variable instead of t (e.g. |2^ 
u = = v/ad in respect with this time variable. Then eq. 



and introduce a comoving velocity 
has the form 



du 
da 



+ (mV, 



>u 



[aa] 



{3a^d)u — Vx(, 



(4) 



Let $ be velocity potential so m = Va;$. We also observe that the combination A = {^H^a^)~^ = 
— (3aa^)~^ does not depend on time in the matter-dominated Einstein-de Sitter Universe. Then 
from eq. (H) we find a general relation between the velocity potential and the Newtonian gravita- 
tional potential ^5 

+ -(V.$)' =-- 



2a (* + ^^: 



(5) 



da ' 2 

This equation is often reffered to as the Bernoully equation. Substituting density p from the right 
hand side of eq. (H) into eq. (|l|), we find the second equation linking scalar potentials and $ 



v. 



0. 



(6) 
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The structure of this equation is ■ S = 0, where the vector S(a;) is the expression in the square 
brackets in eq. (BI). Using the Helmholtz theorem, we find 



An J \x — x 

After some vector algebra, we end up with the formula 



-v.. X S(f' ^ 



(7) 



d_ 
da 



a \ f d^x' 



Thus we obtain two equations (|]), (|]) for two fields $ and 
construct the traceless symmetric tensors a/y = V.r.V^^.^ 



In the Newtonian theory we can 



|5.,V^$andE,^ 



Tensor Ef^ is the gravitational tidal field, and corresponds to the Newtonian limit of the electric 
part of the Weyl tensor. 

Taking derivatives V Xi^Xj from eq. (^, we can get the equation for the Lagrangian time derivative 
along the trajectory -^cr^ plus other local terms only. Taking derivative V^^ of eq. (|^), we obtain 
equation containing the term -^E^ plus other local terms in the left hand side, and non-local 



integral term V^^ .V^-feVa;^ 



/#|.(V..$)(V^- 



-x'\ 



in the right hand side. Apparently, this non-local 
the local terms. It prevents us from obtaining a 



term cannot be reduced to the combination o 
closed set of Lagrangian equations in the Newtonian theory. However, in general relativity there is 
a closed set of Lagrangian equations for the Weyl tensor and (Xjj [jl^ . In general relativity 00 com- 



ponent of the Einstein Eqs. (which contains the Poisson equation as the Newtonian limit) plays 
the role of constraint equation, which is authomatically resolved for arbitrary time t > to if it is 
resolved once at the initial hypersurface t = to and evolution equations are satisfied. Recently this 
fact drew much attention in connection with possible application to the LSS dynamics beyond the 
linear theory |T^, However its implication for the Newtonian theory is unclear, because 
in the Newtonian limit the Poisson equation has to be resolved for each moment t. 
In this contribution I will consider the general method of following the time evolution of PDFs for 
the dynamical systems which obey the closed set of Lagrangian equations. In the case of New- 
tonian gravity, we will use an approximation where basic equations are truncated and might be 
reduced to the closed set of equations containing full time derivatives only. 



3 Truncated Zel'dovich Approximation 

The dynamics of the system which obeys the multiple stream generalization of the basic equations is 
complicated and requires the N-body simulations. For interesting cosmological models such as the 
CDM scenario, the structure formation looks like complicated hierarchical pancaking and clustering 
from very small to large cosmological scales. However the gravitational clustering at sufficiently 
large scales R can be considered in the quasilinear theory in a single stream regime ignoring small 
scale details. For this goal let us apply the Zel'dovich approximation for the smoothed initial 
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gravitational potential filtered by the window function W{R) with filtering scale R: 



oo 

'3, 



0(f;i?) = — / d'(P{x')W{\x-x'\;R). (9) 

ZTT^ JO 

We will call this approach the truncated Zel'dovich approximation, which was used for different 
purposes in papers g, |2|, [gS], pf, [H, [|| 



The mean space separation between galaxies is about Rq ~ 5h~^Mpc. To generate the continuous 
field from galaxies , one has to smooth the survey with the filter exceeding the minimal scale 
Rq. At these scales the density contrast 6p/p less than unity and quasilinear truncated Zel'dovich 
approximation can be applied. The approximated Zel'dovich solution |^ of the basic equations 
describes the gradient mapping of the Lagrangian particles space q into physical Eulerian space x. 
The velocity field is defined by the gradient of the initial velocity potential ^{q). In this approxi- 
mation the follow of particles density field evolves as p{x) = p{q) |||| . From mathematical point 
of view the truncated Zel'dovich approximation means just the Zel'dovich approximation which 
is applied to the truncated initial potential (|) to ensure being in the single stream quasilinear 
regime. In this Section the approximation is formulated in terms of an equation for the gravita- 
tional potential and its derivatives corresponding to the tidal force. 

In the Zel'dovich approximation = —A~^<^, and then the dynamical equation for the velocity 
potential is reduced to the "shortened" eq. (^) without the right hand side: 

^4(V.«)^ = 0. (10) 

We introduce the tensor of the velocity derivatives Sij{x, t) = V x^Vj in the Eulerian space; for the 
potential motion it is reduced to Sij = Vx.'Vxj^- Let A, be its eigenvalues. The field of the Sij{t)- 
tensor evolves in time, its initial value (in the Lagrangian space) coincides with the Lagrangian 
deformation tensor Dij = Vg, V^^ ^o- We will call the S'jj-tensor as the Eulerian deformation tensor. 
It is convenient to use three invariants of the S'jj-tensor: 

^1 = Ai + A2 + A3; J2 = A1A2 + A1A3 -|- A2A3; J3 = A1A2A3. (11) 

Taking the derivatives of $ in respect with V^, we get from eq. (p!0| ) the following set of equations. 
First equation is 

Du du , , ^ , , 

_ . _ + („vj« = 0, (12) 

which is the truncated Euler equation (cf. eq. ^)) describing the free streaming of particles in a 
comoving coordinates ( is the lagrangian derivative along the trajectory). Further derivatives 
of this equation give us the dynamical equations for the invariants of the deformation tensor field 

DJi . ,2 



Da 

DJ2 
Da 

Da 



Jf - 2J2 = 0, (13) 

+ JiJ2-3J3 = 0, (14) 
J1J2 = 0. (15) 
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Also, we can rewrite the continuity equation in new variables. Let q = a^p be the comoving 
density, then from eq. (|ip we have 

^ + QJi = 0. (16) 
Da 

The advantage of the set of equations (p!^)-(p!6|) is that it is the closed system of equations which 
describes the dynamics of the cosmological potential density, velocity and deformation tensor fields 
$, g, u, Jl, J2, J3 as function of time a{t) and the Eulerian position x in the Zel'dovich approxima- 
tion. In the next Section, using the system (p!^)-([T6|), we will derive and solve the kinetic equation 
for the joint distribution function of these fields. 

Let us introduce the principal radii of curvature Ri = of the three dimensional hypersurfaces 
of the potential Then Jj can be expressed through i?j-th. The solutions of eq. (p!^)-([T5|) in 
terms of Ri are extremely simple: 

i?, = Roi + a, (17) 

the principal radii of curvature of the hupersurface of the potential in the Zeldovich approximation 
linearly increase with "time" a{t). It is nothing but the Huygens principle of geometrical optics 
(note the optical-mechanical similarity between the ray propagations and particle trajectories in 
the Zel'dovich approximation I^T]). Using (|l3), (P^ we obtain the time evolution of the fields 
from their initial values (Jqi, J02, Jos) along with the characteristic: 

Joi — 2aJo2 + 3a^ Jo3 , . 

1 - aJoi + aVo2 - a-^Jos 

J02 — 3aJo3 , , 

1 — a J 01 + a^Jo2 — CL Jo3 

= -^ r^r%7 TT^ ^20) 

1 — a J 01 + a J02 — ct Jo3 

^=-^ TTT^T rr^ ^21) 

and $(g, a) = ^o{q), u{q,a) = Uo{q). From the set of eqs. ([T8| ) - (|2l|) we can obtain the reverse 
formulas 

^ Jl + 2a J2 + 3a Vs 
°' l + aJi + aVs-FaVs' ^ ' 

T _ J2 + 3a J3 

•^02 — :r- J — : — 5-7 — : — [^o) 



1 + aJi + a^ J2 -|- a^J; 
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•^03 = . , , , ^\ . , . . , (24) 



1 + aJi + a^ J2 + a'^J; 
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^° ~ 1 + aJi + a2J2 + a3J3' "^^^^ 
Additionally, we can directly get the solution of eq. (0) 

$(a,f) = <l>o(g-) + ^^^. (26) 

The solution (p6D describes the deformations of 3D-hypersurface of the potential $, and is valid 
untill the formation of folds of •I'-hyper surf ace, which corresponds to caustics (pancakes). Taking 
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the derivatives of $ in respect with V^: and Vg, from (^) we get the usual form of the Zel'dovich 
approximation x = q + aVg$o(0- 



4 Kinetic equation for the PDFs 



In the Zel'dovich approximation, we have derived the closed system of the dynamical equa- 
tions (p^ -(|l6|), and found its solutions which describes the evolution of the set of the fields 
$, Q, u, Ji, J2, J3 ( note that, for instance, the set of fields $, g, u does not obey the closed system 
of equations). For the closed system of dynamical equations one can directly derive a so-called 
" kinetic" equation for the joint PDF. In the context of the two dimensional geometric-optics 
problem kinetic equation was derived in [^], p8| . 

Here we derive the one-point joint PDF ^, -u, Ji, J2, J3) of the fields ^, m, Ji, J2, J3 in the 

Zel'dovich approximation. Let us introduce an arbitrary function of these fields /($, £>, "U, Ji, J2, Js)- 
Taking the full time derivative ^ and using the set of equations (p!^)- ([T6|) , we get 

I + ,V.,(S/) - + + (Jl - 2J.)f + - 3J.)I{ + J.J.I^-^ = 0. (27) 

Convolving this expression with W{^, g, u, Ji, J2, J3), and using the fact that the function / is an 
arbitrary one, we obtain the sought kinetic equation for the joint PDF: 

DW d d d d 

-—-JiW- JijriQ^) - TpriJi - 2^2)H^ - jprW^ - SJa)^^ - t7^(Ji JsW^) = 0. (28) 
Ua og oJi 0J2 0J3 

Let Wq{^q, go, Uq, Jqi, J02, J03) be an initial joint PDF of the initial fields ($0) Qo, uq, Jqi, J02, Jos) 
defined in the Lagrangian space. 

Thus, eq. ( ^8]) together with the given initial condition, describes the time evolution of the statistics 
of the field in the problem. Eq. (123) admits a simple analytical solution which can be obtained 



applying the method of characteristics. Using the time evolution (p!8D-(pl|) of the fields along of 
characteristics, we finally obtain the solution of the kinetic equation 

W{^, g, u, Ji, J2, J3; a) = (1 + aJi + 0^2 + aV3)"Vo($o, ^0, ^o, Jou J02, Jos), (29) 



where we have to substitute the expressions (p2D-(p5|) as the arguments of the function Wq. 
Formula (^) is the main result of the paper. It describes the time evolution of the joint PDF of the 
cosmological potential, density, velocity and its derivatives fields in the Zel'dovich approximation, 
from the given (Gaussian or non-Gaussian) initial distribution function Wo- 



Evolution of the PDFs from Gaussian initial fluctua- 
tions 
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5.1 General formalism for Gaussian initial statistics 

In this Section we study the statistics of the continuous cosmological fields evolving from the initial 
Gaussian fiuctuations, which is apparently the most attractive model of primordial perturbation. 
On the other hand, in this case we can make more advanced predictions for the statistics of 
the matter distribution and motion. For the sake of simplicity, we consider the joint PDF of the 
cosmological density, velocity and deformation tensor. For the Gaussian initial conditions we write 
the initial joint PDF as 

Wo{qo,uo, Joi, -^02, ^03) = Qo{uo) ■ S{qo - q) ■ Go{Joi, J02, Jos)- (30) 

The first factor is the Gaussian velocity distribution function. The velocity dependence is factor- 
ized because there is no correlation between Vg$o and other fields involved in (|30| ) for the Gaussian 
fluctuations. The second factor is the initial density distribution function, which corresponds to 
the perfectly homogeneous density distribution qq = g. This is just the formal limit of the Gaus- 
sian density distribution with a — 0. The third factor is the joint distribution function of the 
invariant of the initial deformation tensor. Substituting the form into (^), we get the time 
evolution of the joint PDF from Gaussian initial fiuctuations. 



5.2 Evolution of the velocity PDF from Gaussian initial fluctuations 



The initial PDF of the velocity u{x, h) for the Gaussian field is 



rexp 



2al 



(31) 



where ct„ = OuiR) is the initial dispersion of the velocity m, which depends on the filtering scale 
R. 

From the joint PDF (P^ ) for the Gaussian initial condition ( PU] ) we see that the velocity dependence 
remains to be factorized for an arbitrary moment of time a(t). Integrating the joint PDF over all 
the arguments, except velocity, we get the Eulerian velocity PDF: 



[2vra2]V: 



-exp 



if 



2al 



(32) 



Thus we find that the Eularian velocity PDF Q{u) is time-invariant under the Zel'dovich approx- 
imation, in accordance with paper ||2^, where this result was obtained by a different method, 
and supported by N-body simulations beyond the Zel'dovich approximation. The time-invariance 
means that the Eulerian velocity PDF remains Gaussian with the dispersion defined by the linear 
theory. The distribution of the velocity field u is isotropic. Using eq. (^), we can obtain the PDF 
of the physical velocity v = o-it)^ = adu: it is Gaussian with the dispersion ay{R) defined by the 
linear theory and the smoothing filter. 
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5.3 Evolution of the density PDF from Gaussian initial fluctuations 



To calculate the Eulerian density PDF P{g, a), we have to integrate the joint PDF given by eq. (pOD 
over all of its arguments except density 

P{g, a) = J (fu dJi dJ2 dJz W{g, u, Ji, J2, J3; a) . (33) 

To take the integral (|33|) , we need the joint distribution function of the invariants of the Lagrangian 
deformation tensor Gq^Jqi, J02, Jqs). The joint distribution function of its Lagrangian (initial) 
eigenvalues Aoi is given by 



-^('^Ol, -^02, A03) 



5^/^ 27 
STial 



(Aoi - Ao2)(Aoi - Ao3)(Ao2 - A03) exp 



at 



O 72 



15 



J, 



02 



(34) 



where Jqi, J02 are expressed through Aqj, and (Tj„ = ain{R) is the initial variance of g, which 
depends on the filtering scale R. Then the joint distribution function of the invariants Ji is 



G'( Joi; <^02, J( 



03 J 



5^/^ 27 



1 



0",-, 



15 



o ( 3t/g]^ ^ J{ 



'02 



(35) 



Additionally, we have to integrate ( ^3]) over the allowed region in the ( Jqi, J025 -^03)^space, for 
which all three eigenvalues Aoi are real. To find that region is a non-trivial task, the basic idea is 
outlined in |2^. After that, taking integral over the velocity in eq. (^3]), we get 



P{g,a) 



dJidJ2dJ^ 



[1 + aJi + a?J2 + a^Jz 



]^{Qo — q)Gq{Jqi, J02, Ji. 



03 J 



(36) 



We have to substitute Go from (pSj) and arguments Joi, J02, Jos from (^2D -(P5l) in this formula. 
After some tedious algebra we can reduce the integral ( |36D to the simpler one-dimensional integral 
which has to be performed numerically 



P{g,a) 



AnNsg^a'^ ^3(^)1/3 V / V 



(37) 



Pnis) 



cos 



2/ .X 1 

-In — liTT H — arccos 

3^ ' 3 



5Af 



(38) 



where the only parameter a{t) = a(t)(jj„ is the standard deviation of the density fluctuations g/g 
in the linear theory, and Ng is the mean number of streams, A^^ = 1 in the single stream regime. 
The expression (37) was derived by the different method earlier in ||2^, p9| . 

For cr ^ 1 the expression (|37| ) is reduced to the Gaussian distribution P{g) = (27rcr^)^^/^exp^(| — 1) 

For 0" — >■ we formally get P{g) S{g — g), in accordance with ([30|). The density PDF P{g,a) 
calculated from (|37D is plotted in ||26l, ||29| for different values of parameter a. The matter is 
evacuated from the underdense regions with g < g, resulting in voids which expand in time and 
tend to occupy a larger fraction of volume, as well as formation of anisotropic collapsed dense pan- 
cakes which tend to occupy a smaller fraction of volume. The positive density contrast can reach 
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any large value while the negative density contrast is restricted by p > 0. Hence the probability 
function P{g,a{t)), meaning the fraction of volume with a given value of density, is expected to 
be very non-Gaussian even in the quasilinear stage. At large densities g ^ g PDF from ( pT)) has 
the pancakes induced asymptota oc g~^, which affects such values as skewness, 5*3^'' = 4 in this 
approximation instead of the actual value ~ -y — + 3). However, the final smoothing (or 
the adhesion) regularizes high-density asymptota of P{g), see 0. Formula (pTf ) very well 
describes actual P{g) in the practically interesting region < ^ < several g for cr < 1/2. 



5.4 Evolution of the deformation tensor PDF from Gaussian initial 
fluctuations 



The density and velocity PDFs were considered earlier in the literature. Let us consider the new 
statistics: PDF of the Eulerian deformation tensor. The time evolution of the deformation tensor 
gives us alternative approach to the LSS 0. To calculate its PDF, we have to integrate the joint 
PDF over all the arguments except Jj-th: 



G( Ji, J2, J3; ^) ^ j ^^"^ ^{.9^ Ji^ J2, J3', a) ■ 



(39) 



Carrying out this integral and switching from invariant Jj-th back to Aj-th, we obtain the joint 



distribution of the Eilerian eigenvalues of the deformation tensor 



M(Ai,A2,A3;a) 



55/2 27 



(Ai-A2)(Ai-A3)(A2-A3) 
((l + aAi)(l + aA2)(l + aA3))' 









exp 







, (40) 



where Jqi and J02 must be expressed through Aj via eqs. (p2D -(^5D. In the limit a ^ this distri- 
bution is reduced to the initial Doroshkevich's distribution ( ^4|) based on the Gaussian statistics. 
As cr = a{t)ain is growing, the Eulerian distribution M(Ai, A2, A3; a) rapidly departs from the ini- 
tial Lagrangian distribution, and might be an interesting discriminative test for the initial statistics. 



6 Properties of the density PDF 



6.1 PDF from the Edgeworth perturbation series 



In this Section we consider some properties of the cosmological density PDF using other approaches 
besides the Zel'dovich approximation. In the case of the weakly non-linear dynamics {S < 1) 
when slight departure from the initial Gaussian distribution is expected, one can use the general 
decomposition series around the gaussian PDF ||2^. Inspired by the paper |^ on the weakly 
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non-linearities on statistical distributions in the theory of two-dimensional sea waves, we suggest 
similar decomposition for the cosmological density P{S) (in the three-dimensional case). The 
so-called Edgeworth form of the Cram-Charlier series for density distribution function reads as 



P{6) 



(27ra2)V: 



-exp 



2a2 



1 + cr ■ — ■ _H3 + cr ■ — Ha H H( 

6 V24 72 ^ 



+ 



(41) 



where Hn{^) is the Hermit polynoms. Numerical coefficients Sn can be calculated in the pertur- 
bation series [|I|. We found [0] that a few iterations of the expansion (^) reproduce the peak of 
P((5) in the interval of \5\ < 1/2 around it relatively well for small parameter a < 1/2. It fails to 
reproduce P{5) at larger \5\ because the series (^Tj) is an asymptotic expansion. We understand 
that similar decomposition was independently suggested in [p3 . 



6.2 Mystery of the log-normal distribution 



As it was noted in [jl3l, ||2T 



the lognormal distribution 

(ln^-/i/)2 



(27ra2)i/2 



exp 



2af 



1 

Q 



(42) 



is a good approximation of the actual P{q), fii = In fi — O.bcrf, erf = ln(l + o"^//!^). However, this 
fit was checked for CDM model for moderate a ~ 1/2 only. The question arises, if the log- normal 
distribution is a universal form of P{q) due to the non-linear dynamics of the cosmological system, 
or it is just a convenient fit for particular cosmological models in some intermediate regime? We 
argue p| that the log-normal PDF is not a universal form but is close to the actual P{g) for some 
cosmological models (including CDM) for moderate a. 

One can use the Edgworth seria (^1]) for P{g). From (|4^) we find the skewness of the log-normal 
distribution S^^°^^ = 3 + cr^. On the other hand, the skewness of the filtered density field is 



S: 



_ 34 
3 - - - 

equation 



in. 



+ 3). To be the distributed lognormally, density field has to satisfy the following 



3 -t- cr = y 



in. 



+ 3). 



(43) 



This is a necessary (but not sufficient) condition for the a — riejj-dependence to correspond to 
the perfect log-normal distribution. Clearly, it is not a general condition for the density field in 
arbitrary cosmological models. For instance, for CDM model a rough approximation (for interval 
—2 < Heff < —1) is aijieff) ~ 1.4 ■ {—Ueff — 0.85). This aijieff) dependence for CDM model is 
close to that of log-normal distribution given by eq. (|43|) for moderate a ~ 0.4 — 0.6 only, and 
departs from log-normal formula for small and large a. To support this conclusion, the system- 
atic comparison of P{q) against the log-normal distribution for different a{neff) is needed. Thus, 
some "lognormalish" features in the observed density distribution mean that realistic cosmological 
model is close to the CDM one. 
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